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Abstract 

A nucleation model for the breakdown phenomenon in an initial non-homogeneous 
free traffic flow that occurs at an on-ramp bottleneck is presented. This model is 
in the context of three-phase traffic theory. In this theory, the breakdown phe- 
nomenon is associated with a first-order phase transition from the "free flow" phase 
to the "synchronized flow" phase. In contrast with many other nucleation models for 
phase transitions in different system of statistical physics in which random precluster 
emergence from fluctuations in an initial homogeneous system foregoes subsequent 
cluster evolution towards a critical cluster (critical nuclei), random precluster oc- 
currence in free flow at the bottleneck is not necessary for traffic breakdown. In 
the model, the breakdown phenomenon can also occur if there were no fluctuations 
in free flow. This is because there is a permanent and motionless non-homogeneity 
that can be considered a deterministic vehicle cluster localized in a neighborhood 
of the bottleneck. The presented nucleation model and a nucleation rate of traffic 
breakdown that follows from the model exhibit qualitatively different features in 
comparison with previous results. In the nucleation model, traffic breakdown nucle- 
ation occurs through a random increase in vehicle number within the deterministic 
vehicle cluster, if the amplitude of the resulting random vehicle cluster exceeds some 
critical amplitude. The mean time delay and the associated nucleation rate of traf- 
fic breakdown at the bottleneck are found and investigated. The nucleation rate of 
traffic breakdown as a function of the flow rates to the on-ramp and upstream of the 
bottleneck is studied. Boundaries for traffic breakdown in the diagram of congested 
patterns at the bottleneck are found. These boundaries are qualitatively correlated 
with numerical results of simulation of microscopic traffic flow models in the context 
of three-phase traffic theory. 
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1 Introduction 



Empirical observations of freeway traffic made in various countries show that 
the onset of congestion in an initial free flow is associated with an abrupt 
decrease in vehicle speed. This traffic breakdown called the "breakdown phe- 
nomenon" occurs mostly at freeway bottlenecks, in particular on-ramp bottle- 
necks. The traffic breakdown is accompanied by a hysteresis effect (see refer- 
ences in the reviews [1,2], the book [3], and the conference proceedings [4,5,6]). 
The breakdown phenomenon has a probabilistic nature [7,8,9]: At the same 
on-ramp bottleneck, traffic breakdown is observed at different flow rates in 
different realizations (days). The probability of the breakdown is a strong 
increasing function of flow rate downstream of the bottleneck [8,9]. 

Most microscopic, macroscopic, probabilistic, and other models of freeway 
traffic explain the onset of congestion in free flow by moving jam emergence 
(see, e.g. [10,11,12,13,14,15] and references in the reviews [16,17,18,19], as 
well as the conference proceedings [20,21,22,23,24]). In particular, in mod- 
els of an on-ramp bottleneck moving jam(s) occurs spontaneously in free 
flow at the bottleneck when the flow rate upstream of the bottleneck is high 
enough and the flow rate to the on-ramp increases gradually beginning from 
zero [17,18,25,26,28,27]. However, the fundamental model result that the onset 
of congestion in free flow on a homogeneous road and at freeway bottlenecks is 
associated with spontaneous moving jam emergence [10,13,14,15,25,26,28,27,16,17,18,19] 
is in a serious conflict with empirical evidence [29,30,3]. 

Consequently, in 1996-1999 Kerner introduced three-phase traffic theory (see [3] 
for a review). In this theory, there are three traffic phases: free flow, synchro- 
nized flow, and wide moving jams. In accordance with empirical investigations 
of phase transitions, in this theory moving jams do not emerge spontaneously 
in free flow. Rather than moving jam emergence, a phase transition from free 
flow to synchronized flow (F— >S transition for short) governs the onset of con- 
gestion in free flow [29,30]. A first-order F— »S transition postulated in three- 
phase traffic theory [29] discloses the nature of the breakdown phenomenon at 
freeway bottlenecks found in empirical observations [1,2]. In other words, the 
terms "F-^S transition", "breakdown phenomenon" , "speed breakdown" , and 
"traffic breakdown" are synonyms. The first microscopic models in the con- 
text of three-phase traffic theory are stochastic models [31,32]. These models 
exhibit phase transitions as well as all types of congested patterns found in 
empirical observations [31,32,33,34,3]. Recently, some new microscopic models 
based on three-phase traffic theory have been developed [35,36,37], which can 
show some congested pattern features found earlier in [31,32]. 

One of the important methods for a study of phase transitions in non-linear 
distributed multiple-particle systems is a probabilistic theory based on an 
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analysis of a master equation (e.g., [38,39,40]). First probabilistic theories for 
traffic flow based on a master equation for a random vehicle cluster have 
been introduced by Mahnke et al. [13,14] and further developed by Kiihne et 
al. [15,41] (see the recent review by Mahnke et al. [40]). As usual for other 
metastable systems of statistical physics, for an initial metastable traffic flow 
a well-known two-well potential nucleation problem arises from the master 
equation, which analysis is made based on general well-known methods of sta- 
tistical physics [38,39]. One of the main results of this analysis is the nucleation 
rate for the critical vehicle cluster (critical nuclei) whose occurence leads to a 
phase transition within the initial metastable state of traffic flow. As usual for 
other metastable systems of statistical physics [38,39], in the metastable traffic 
flow the nucleation rate for a phase transition is an exponential function of a 
control parameter, flow rate (or vehicle density) for traffic flow [13,14,15,41,40]. 
Rather than these common well-known results, an interest for the field has a 
nucleation model for a metastable traffic flow. Both the model and associated 
dependencies of the nucleation rates for phase transitions on traffic variables 
and control parameters should be adequate with real traffic flow features. 

In [13,14,15], models for vehicle cluster nucleation in an initial spatially ho- 
mogeneous traffic flow on a homogeneous circular road are introduced (see, 
however, Sect. 4.2). One of the basic assumptions of these models is that in 
an initial homogeneous free flow firstly random precluster should emerge from 
fluctuations. In other words, this precluster foregoes subsequent cluster evolu- 
tion towards a critical cluster (critical nuclei) whose growth leads to a phase 
transition. The rate of precluster emergence w + (0) in traffic, in which initially 
no vehicle cluster exists, can be different from the attachment rate of cluster 
evolution, when a random cluster with n > 1 vehicles already exists at the 
road [13,14,15,40]. 

A first nucleation model based on the master equation for traffic breakdown 
at an on-ramp bottleneck, i.e., in an open traffic system has been suggested 
by Kiihne et al. [41] and Mahnke et al. (Chap. 17 of Ref. [40]). As in the 
models of homogeneous road [13,14,15], in this model at given flow rates to 
the on-ramp q on and on the main road upstream of the on-ramp q in random 
vehicle cluster occurrence and evolution that can lead to traffic breakdown 
are realized only after a vehicle precluster consisting of one vehicle randomly 
appears at the bottleneck. The rate of precluster emergence is equal to the 
flow rate to the on-ramp [41,40]: w + (0) = q on . This idea about foregoing 
precluster emergence is associated with a basic model assumption that at q on = 
no vehicle cluster can randomly appear, specifically no traffic breakdown is 
possible [41,40]. Probably, this basic model assumption leads to the nucleation 
rate for traffic breakdown that is proportional to q OI1 [41,40]. However, in real 
traffic flow both flow rates q on and q in exhibit random fluctuations that can 
cause cluster emergence. 
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Moreover, whereas the basic assumption about the necessity of precluster 
emergence is physically justified for a homogeneous road [13,14,15], this is 
not the case for for traffic breakdown at the bottleneck. This is because in 
contrast with the model of a homogeneous road [13,14,15], an initial state of 
free flow at the bottleneck at q OIl > and q m > is non-homogeneous re- 
gardless of fluctuations [3]. This means that even there were no fluctuations 
in free flow at the bottleneck, nevertheless free flow is non- homogeneous in a 
neighborhood of the bottleneck. This is because two different flows permanent 
merge within the merging region of the on-ramp - the flow with the rate q on 
and the flow with the rate q in . 

It has been shown in three-phase traffic theory [3] that due to the merging of 
these flows a steady state of free flow at the bottleneck is non-homogeneous: 
A permanent and motionless local perturbation (deterministic perturbation) 
occurs at the bottleneck. Within this perturbation that can be considered a de- 
terministic vehicle cluster the speed is lower and density is greater than down- 
stream of the cluster. This deterministic cluster can lead to traffic breakdown 
even if there were no random fluctuations at the bottleneck [42,3]. However, in 
the nucleation model [41,40] the master equation is formulated for probability 
of a vehicle cluster, which can occur due to fluctuations only. We see that 
for an initially non-homogeneous traffic flow, which occurs at the bottleneck, 
another physical approach should be applied. This is due to at least two rea- 
sons mentioned above: (i) Regardless of fluctuations, there is already a vehicle 
cluster at the bottleneck that growth can lead to traffic breakdown: No vehicle 
precluster emerging from fluctuations is necessary for traffic breakdown, (ii) 
When influence of fluctuations on traffic breakdown is considered, fluctuations 
caused by both flow rates q on and q m should be taken into account. 

Thus, in contrast with the nucleation model at the bottleneck of Ref. [41,40] 
a nucleation model that is adequate with empirical results for traffic break- 
down at the bottleneck should consider vehicle cluster evolution in an initially 
non- homogeneous free traffic. This random vehicle cluster should include both 
the total vehicle number within an initial deterministic cluster and addition 
vehicles within the cluster, which occur due to fluctuations. 

However, there are no nucleation models and associated probabilistic the- 
ories for an initially spatially non-homogeneous free flow. Moreover, struc- 
ture non-homogeneities play also a very important role for phase transitions 
in many other multiple-particle systems of statistical physics like granular 
flow, biological physics, reaction-diffusion systems, etc. (see references, e.g. 
in [20,21,22,23,24,43,44]). Thus, it seems important to derive a probabilistic 
theory for traffic breakdown in non-homogeneous traffic flow. 

In this paper, a nucleation model for the probabilistic breakdown phenomenon 
in an initially spatially non- homogeneous traffic flow at an on-ramp bottleneck 
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is presented. The article is organized as follows. In Sect. 2, the nucleation model 
is considered. Nucleation rate and the mean time delay for traffic breakdown 
that result from this model are studied in Sect. 3. Results of the paper as well 
as their comparison with earlier nucleation models for traffic breakdown of 
Ref. [13,14,15,41,40] are discussed in Sect. 4. 



2 Nucleation model of traffic breakdown in initially non-homogeneous 
free flow at on-ramp bottleneck 

2. 1 Non-homogeniety in free flow at bottleneck as reason for vehicle cluster 

In accordance with three-phase traffic theory [42,3], in a nucleation model we 
assume that the breakdown phenomenon at an on-ramp bottleneck is asso- 
ciated with an initial non-homogeneity of free flow at the bottleneck. This 
non-homogeneity exists at the bottleneck even in a hypothetical case in which 
there were no fluctuations in traffic flow. In this hypothetical case, this non- 
homogeneity can be considered a deterministic (permanent) local perturba- 
tion localized at the bottleneck. The non-homogeneity of an initial free flow 
at the bottleneck is caused by two flows, which merge at the bottleneck: (i) 
An on-ramp inflow with the rate q on . (ii) A flow on the main road upstream 
of the bottleneck with the rate q in . This flow merging occurs permanent and 
on the same freeway location (within an on-ramp merging region). For this 
reason, the non-homogeneity of free flow at the bottleneck is motionless and 
permanent (Fig. 1 (a)). 

To explain features of this non-homogeneity, note that at a given high enough 
flow rate q in in free flow on the main road upstream of the bottleneck, vehicles 
that merge from the on-ramp onto the main road force the vehicles on the main 
road to decelerate in the vicinity of an on-ramp merging region. This deceler- 
ation leads to a local decrease in speed and consequently to a local increase 
in density in the vicinity of the bottleneck. As a result, a local perturbation, 
i.e., non-homogeneity in free flow appears. 

If no fluctuations in free flow are considered, as mentioned above the non- 
homogeneity in free flow in the form of a deterministic local perturbation 
occurs at the bottleneck. The speed v^fj e and density p\fj e within this deter- 
ministic perturbation correspond to the conditions 



where -y( free ) and p( frce ) are the average vehicle speed and density in homo- 
geneous free flow on the main road downstream of the perturbation (Fig. 1 
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Fig. 1. Explanation of the basis of nucleation model: (a) Sketch of an on-ramp bot- 
tleneck, (b, c) - Z-shaped speed-flow (a) and the associated S-shaped density-flow 
characteristics for an F— >S transition, (d, e) - Simplified Z-shaped speed-flow (d) 
and S-shaped density-flow characteristics (e) related to (b) and (c), respectively. 

(a)). Because the deterministic local perturbation is permanent and motion- 
less, at given g in and q on , the total flow rate does not depend on the spatial 
co-ordinate, i.e., the speed and density on the main road satisfy the following 
condition: 



n — ? , (free) (free) 



V (B) (B) 
17 freer free' 



(2) 



where 



Qsum — Q'm 4" Qon- (3) 

The inhomogeneous steady state within the deterministic local perturbation 
can be considered "deterministic vehicle cluster" in free flow localized at the 
bottleneck or "deterministic cluster" for short. 

When q on is a given value and the total flow rate q sum increases, then the vehicle 
speed VffJ e within the deterministic cluster decreases and in accordance with 
(2) the associated density p^J e increases. However, this increase is limited by 
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some critical density p^fj e = preterm, fs (Fig- 1 ( c )) within the deterministic 
cluster associated with a critical flow rate 



Q-sum *3determ, FS ' 

(4) 



After this critical deterministic perturbation is reached, the further increase in 
<?sum leads to deterministic traffic breakdown at the bottleneck causing sponta- 
neous synchronized flow emergence at the bottleneck. The critical determin- 
istic perturbation can be considered "critical deterministic vehicle cluster" 
("critical deterministic cluster" for short) or "deterministic nuclei for traffic 
breakdown" . After the critical deterministic cluster is reached, deterministic 
traffic breakdown occurs at the bottleneck even if there were no random per- 
turbations in traffic flow at the bottleneck. 

Random perturbations within the initial deterministic cluster can cause ran- 
dom traffic breakdown (F^S transition) at the flow rate 

Qsum Qdeterm, FS ' (5) 



i.e., before the deterministic nuclei for traffic breakdown is reached. In this 
case, random traffic breakdown nucleation can occur at the bottleneck (arrows 
labeled F— >S in Fig. 1 (b, c)). This is realized, if through a random increase in 
vehicle number within the initial deterministic cluster, the amplitude of the 
resulting "random vehicle cluster" ("random cluster" for short) exceeds some 
critical amplitude associated with a critical density within the random cluster 
Pct^fs (Fig- 1 (c)). The random cluster with the critical density Pcf^s can 
be considered "critical random vehicle cluster" at the bottleneck ("critical 
random cluster" for short) or "random nuclei for traffic breakdown". If the 
amplitude of a random cluster is smaller than the critical one, the random 
cluster decays towards the initial deterministic cluster. 



2.2 Master equation 



2.2.1 Basic assumptions for master equation 

There are the following basic assumptions for the nucleation model in an initial 
non-homogeneous free flow at the bottleneck: 

(i) There is no vehicle precluster, which random occurrence through a random 
decrease in speed of one of the vehicles in the initial free flow should be nec- 
essary for vehicle cluster emergence and subsequent cluster evolution. Traffic 
breakdown occurs within a deterministic vehicle cluster that is motionless and 
permanently exists in a neighborhood of the bottleneck. 
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cluster size growth dissolution 

Fig. 2. Schematic illustration of cluster transformation in an initially 
non-homogeneous free flow at the bottlneck. 

(ii) A master equation should describe probability p for random vehicle cluster 
evolution in which the cluster size N, i.e., the total vehicle number within the 
motionless vehicle cluster randomly changes due to fluctuations in a neighbor- 
hood of the deterministic cluster (Fig. 2). The size of the deterministic cluster 
jy(dctcrm) does not depend on fluctuations in traffic flow. 

(hi) The attachment rate w + onto this vehicle cluster is not a function of the 
cluster size N, i.e. w+(0) = w + (N) = q in + q on . 

(iv) A dependence of the detachment rate u>_ from a cluster is a non-linear 
function on N, which consists of at least two branches associated with the de- 
terministic cluster and a critical vehicle cluster (nuclei for traffic breakdown). 

(v) The on-ramp inflow and the flow upstream of the bottleneck make a differ- 
ent influence on the cluster size. To take into account this assymetric behavior 
in the model, the detachment rate W- also depends on q on . 

The assumption (i) can be explained by existence of a deterministic vehicle 
cluster at the bottleneck. This cluster can lead to traffic breakdown at the 
bottleneck even there were no fluctuations in traffic flow, i.e., no precluster 
is necessary for vehicle cluster occurrence in an initially steady state of non- 
homogeneous free flow at the bottleneck. To explain the assumption (ii), note 
that real random fluctuations can cause either a decrease or an increase in 
the cluster size N and, respectively, either a decrease or an increase in vehi- 
cle density within the cluster. Real random fluctuations, which decrease the 
cluster size N, are associated with an increase in speed within the cluster 
(Sect. 2.1). Therefore, these fluctuations can prevent the deterministic break- 
down. In contrast, random fluctuations, which increase the cluster size N, i.e., 
decrease the speed within the cluster, can cause traffic breakdown before the 
deterministic nuclei for traffic breakdown is reached, i.e., when the condition 
(5) is satisfied. Regardless of the cluster size N, the attachment rate into the 
cluster is determined by fluctuations in both flow rates g; n and q on . This ex- 
plains the assumption (iii). Grounds for the assumptions (iv) and (v) appear 
in Sect. 2.2.2. 
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We consider the dynamics of the total vehicle number N within a vehicle 
cluster localized at the bottleneck (dashed region in Fig. 1 (a)). It is assumed 
that within the cluster free flow is spatially non- homogeneous, specifically the 
region of spatial cluster localization is bounded upstream and downstream 
with homogeneous free flows. The total number of vehicles N within the cluster 
can either increase or decrease over time randomly in comparison with a value 
N = jy( determ ) for the case in which the deterministic cluster exists at the 
bottleneck only. 

In accordance with the above basic assumptions, the probability p(N, t) to 
find N vehicles within the cluster at the bottleneck reads as follows 



= w+(N - l)p(N -l,t)+W-(N+ l)p(N + 1, t) 

~[w+(N) + w_(N)]p(N, t), at N > 0, (6) 



^M = w _(i)p(i >t )_ u , + ( )p(0,0, atiV = 0, (7) 
with the boundary condition 



w-(0) = 0, (8) 

In accordance with the assumption (3), the vehicle attachment rate w + is 
independent of N, i.e., 



W+ = q S um- (9) 



2.2.2 Vehicle detachment rate from cluster 

The vehicle detachment rate W- is obviously equal to the outflow rate from 
the cluster 



v>-(N) = $!? ) m (io) 

To understant a qualitative shape of the function q^^ e \N), note that in 
accordance with three-phase traffic theory [3], speed states within the de- 
terministic cluster VffJ e (q sum ), the speed within the critical random cluster 
v cr Fsfeum), along with a 2D synchronized flow speed states [3], together form 
a Z-shaped speed-flow characteristic for traffic breakdown (Fig. 1 (b)). The 
associated density-flow characteristic, which consists of density states within 
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the deterministic cluster p^teum) , the density within the critical random 
cluster Pcf ^Fsl^sum)) along with a 2D synchronized flow states, has obviously 
a S-shaped form (Fig. 1(c)) 1 Due to an F^S transition, fast cluster growth 
at the bottleneck occurs leading to congested pattern formation, i.e., either a 
synchronized flow pattern (SP) or a general pattern (GP) appears upstream of 
the bottleneck [3]. However, the nucleation effect leading to traffic breakdown 
and its characteristics are fully independent of possible congested patterns re- 
sulting from this F^S transition. For this reason, we can average the infinity 
of synchronized flow states (dashed regions in (b, c)) for each given flow rate 
g sum to one speed ug^ avcr (g SU m) and to one density pg, avcr (g S um) (Fig. 1 (d, 
e)) 2 . A consideration of the resulting congested patterns is beyond the scope 
of this article. 

In the model, it is assumed that the shape of the characteristic q^^ e \N) 
(Fig. 3) follows from the S-shaped density-flow characteristic of three-phase 
traffic theory (Fig. 1 (e)): The characteristic ?down^(-^0 has at least two differ- 
ent branches q^°\N) labeled 7V( determ ) an d N c in Fig. 3 (a). These branches 
are related to the vehicle number ranges, respectively, given by the conditions 

< N < N d (11) 

and 

N d <N< N s . (12) 

The branches AK determ ) and N c in Fig. 3 (a) are associated with the density 
branches p[f] e and p^ps °f the S-shaped density-flow characteristic in Fig. 1 
(e), respectively. The branch ]\[( dctCTin ) [ s associated with the case in which at a 
high enough flow rate q Bum and the on-ramp flow rate q on > the deterministic 
cluster exists at the bottleneck. The branch N c is associated with the case in 
which the critical random cluster whose growth leads to an F^S transition 
occurs at the bottleneck. 



Within the deterministic cluster the speed v^' and density p^' shown in Fig. 1 
(b-e) are connected by the formula (2). The number of vehicles within the cluster 
N = f pdx, where integration is performed over the region of cluster localization. 
2 Concerning synchronized flow states in the vicinity of the bottleneck associated 
with a resulting congested pattern, we should note that assumptions used in our 
nucleation model are satisfied only for a localized SP (LSP) [3]. An LSP is an SP 
whose downstream front is fixed at the bottleneck. The upstream front of the LSP 
is localized on the main road at some finite distance upstream of the bottleneck, i.e., 
the width (in the longitudinal direction) of the LSP is always limited. Synchronized 
flow states at the bottleneck are drawn in Fig. 1 (b, c) and in other illustrations 
only for the case in which the congested pattern is an LSP. 
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Mdown (a) 




N, N d N 2 N S N 3 N 




Fig. 3. Qualitative dependencies of the outflow rate q down on the total vehicle 
number N within the cluster localized at the bottleneck (a), and possible depen- 
dencies of the N-shaped function (N) on q on for two different values q on (b); 

if (1) of (2) ^ (1) 

curve 1 for q on = , curve 2 for q on = <?V > q 6 n ■ 

In addition, from the S-shaped density-flow characteristic (Fig. 1 (e)) can be 
seen that for the case in which an LSP results from an F— >S transition, at 



N > N s , (13) 

there is a third branch 7V (syn ' ) on the characteristic (JV) (Fig- 3 (a)) 

associated with the branch aver for averaged synchronized flow states in 

Fig. 1 (e). In this case, q^n^ (N) (10) is a N-shaped flow-vehicle-number 
characteristic. 

At the critical point N = N d at which the branches ]\[( determ ) and N c merges, 
the function q d h ^^ (N) has its maximum point. At the threshold point N = N s 
at which the branches iV c and N^ syn ^ merges, the function q 1 ^^ (N) has its 
minimum point. 

Quantitative characteristics of the N-shaped function q d b ^^ c \N) (e.g., values 
N d and N s ) can depend on the flow rate to the on-ramp q on (assumption (v) 
in Sect. 2.2.2). This is because the on-ramp inflow and the flow upstream of 
the bottleneck can make a considerable different influence on the cluster size 
and the outflow rate from the cluster ■ m particlular, it can turn out 
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that at the same N the greater o on , the more difficult for vehicles to escape 
from the cluster, i.e., the less qj^n^ This is confirmed by microscopic 
simulations [32] and reflected in Fig. 3 (b) in which it is assumed that the 
greater q on , the greater N d and N s are. Thus, in a general case instead of (10) 
we should use 

™-(N) = q £Z le) (N, q on ). (14) 



A possible impact of the flow rate q on on quantitative characteristics of the 
mean time delay for an F^S transition is discussed in Sect. 4. 

Through the use of the basic assumptions (i)-(v) of Sect. 2.2.1 and the chosen 
shape of the function (10) discussed above (Fig. 3 (a)), critical cluster occur- 
rence describes an F— >S transition at the bottleneck. There are two reasons for 
this statement: (i) A vehicle cluster is motionless, i.e., fixed at the bottleneck. 
This is related to the definition of synchronized flow whose downstream front 
is usually fixed at the bottleneck, whereas the downstream from of a wide 
moving jam propagates trough any bottleneck while maintaining the mean 
downstream jam velocity [3]. (ii) The shape of the chosen function (10) in the 
nucleation model associated with this motionless cluster (Fig. 3 (a)) follows 
from a Z-shaped speed-flow characteristic for traffic breakdown (Fig. 1 (b)) for 
an F— »S transition at the bottleneck found in a microscopic traffic flow theory. 
In this theory has been shown that if these two requirements are satisfied, then 
rather than an F-^J transition (moving jam emergence) an F— *S transition, 
i.e., synchronized flow emergence occurs at the bottleneck. As follows from 
this microscopic theory, after the F-^S transition has already occurred, mov- 
ing jams can emerge in this synchronized flow. However, the nucleation model 
describes traffic breakdown, i.e., an F-^S transition, specifically the rate of 
traffic breakdown (synchronized flow) nucleation only. 

In addition, it should be noted that the branch for synchronized flow vj^ aver in 

Figs. 1 (b-e) and the associated branch A( syn ) on the characteristic ^down^ (N) 
(Fig. 3 (a)) follow from the microscopic traffic theory, rather than from the nu- 
cleation model. This branch, which is shown only with the aim of a qualitative 
illustration of a possible traffic flow state after synchronized flow nucleation, 
has no influence on the nucleation rate of an F— >S transition. 



2.3 Steady states 



Steady states of vehicle number N at given q m and q on are associated with 
solutions of the equation 
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w + = w-(N). 



(15) 



In accordance with (9), (14), they are found from the condition 



gsum = ?down e) (^ Qon)- ( 16 ) 

As can be seen from Fig. 3 (a), at given flow rates q on and q\ n that satisfy the 
condition 

<?th < Qsum < <?determ, FS 

(17) 



there can be at least two steady states: N — Aq associated with the de- 
terministic cluster and N = N 2 associated with the critical random cluster. 
These steady states are the roots of Eq. (16), i.e., they are associated with 
the intersection points of the horizontal line q^n^ = Qsum with the branches 
jy(dctcrm) an( j ^ Q j? characteristic ^own^ (N, q on ) (Fig. 3 (a)), respectively. 
In addition, if an LSP occurs as a result of an F— >S transition, then there is 
a third root of Eq. (16), N = A3, associated with the intersection point of 
the horizontal line ^down'^ = q S um with the branch jV (syn ' ) of the characteristic 

If the flow rate g sum increases, then the critical vehicle number difference within 
the cluster 

AJV C = N 2 -N 1 (18) 



decreases. This critical vehicle number difference is associated with the vehicle 
number difference within the critical random cluster and within the initial 
deterministic cluster at the bottleneck. The growth of the critical random 
cluster leads to traffic breakdown at the bottleneck. 

At the critical flow rate (4), we get AiV c = 0: The steady states Ni and 
N 2 merge into one point with the critical vehicle number N — at which 
firm, fs = 9down Ie) (^d, <?on). At g sum > q^ erm> Fs the deterministic traffic 
breakdown should occur even if there is no random increase in the vehicle 
number within the initial deterministic cluster at the bottleneck. 

If the flow rate g sum decreases gradually, then the threshold flow rate 

<?sum = qS ] (19) 



is reached at which the steady states N 2 and A3 merge into one threshold 
steady state N = N s at which q$ = q ( ^ e) {N s , q on ). 
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Fig. 4. Qualitative shape of the potential <&(N) (22) for different flow rates g sum : 
curves 1,2, and 3 are related to the corresponding flow rates qium, qium, and qi^ m 



(3) ! ! 

satisfying the condition Qsum ^ <?sum 
3 Nucleation rate of traffic breakdown at bottleneck 



As follows from the analysis of the model (6)-(9), (14) (Appendix A), in the 
flow rate range (17) the mean time delay of an F— »S transition at the bottle- 
neck is 



Tpg' mean) = Cexp{A$}, 
where a potential barrier 



(20) 



A$ = $(iV 2 ) - $(iVi), 
the potential $(N) is 



(21) 







W— (n) 

Ul_|_ 



at JV > 0, 
at iV = 0, 



(22) 



C = 2n(w'_(N 1 ) | w'_(N 2 ) 



(23) 



w'_(N) = dw-/dN. Respectively, the nucleation rate for traffic breakdown at 
the bottleneck is 



G 



(B) = 1 

FS rp(B, mean) 



C^exp-f- A$}. 



(24) 



To find a qualitative shape $(iV) (22) (Fig. 4), a change in <&(N) between two 
neighboring points N and N — 1 that equals 



5$(N) = $(N) - <$>(N - 1) = In 



W-(N) 



w+ 



(25) 
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can be used. The value 5<&(N) (25) becomes zero at the maximum and min- 
imum points of the function <&(N), i.e., at the roots of Eq. (15) that are the 
points N = N u i = 1,2,3 discussed above (Fig. 3 (a)). The value 5$(N) > 
at w-(N) > w + , i.e., at points of the curve w-(N) above the horizontal line 
q = <?sum in (Fig. 3 (a)). In contrast, 5<&(N) < at w-(N) < w + , i.e., at points 
of the curve W-(N) below the horizontal line q = q sum - 

It can be seen from (20) that the mean time delay for traffic breakdown de- 
creases exponentionally with increase in potential barrier A<£> (21). If in Fig. 4 
the total flow rate increases from q sum = q^ m to q sum = q^ m , which is close to 
the critical flow rate (4) for deterministic traffic breakdown, then the potential 
barrier A$ (21) decreases from A<3>i to A$2- 

In contrast, if the total flow rate decreases from q sum = q^ m to g sum = qj$ m , 
which is close to the threshold flow rate g th (19) for random traffic breakdown, 
then the potential barrier A$ (21) increases from A$i to A<3> 3 (Fig. 4). At 
the threshold point g sum = q^ (19), the potential barrier A<l>(iV) reaches the 
maximum value 



A$ = $(Ag - $(7V th ), (26) 

where N th = Ni at q smn = q^ . As a result, the mean time delay Tpg' mean ' ) 
(20) strongly increases as g sum approaches the threshold point . Under the 
condition 

Qsum < qS ] (27) 



no traffic breakdown at the bottleneck regardless of a random increase in the 
vehicle number within the cluster is possible at the bottleneck. 

If in the vicinity of the critical vehicle number the function w_(N) (14) can 
be approximated by a parabolic function of N, then the following approximate 
formula can be derived from (20) (Appendix A): 

^(B, »„„) V2*N d ( 1 + Ay 2 \ iVvTm , 4V2jV d A;/^ 

FS "^, re « i A )./ 2 l 1 _A^J ex n- VU )' 



where 

= -(N 2 d 2 lnw./dN 2 )\ N=Nd (29) 
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is a dimensionless value of the order of 1, 

(B) 

a _ ^determ, FS Qsura , , 

^deterrn, FS 



i.e., A c is the relative difference between the critical flow rate g determi F s f° r 
the deterministic F— >S transition and the total flow rate g sum (3). If in (28) 
A c <C 1, then we get 

7$ — > = m ^ exp f^T) . (31) 



Respectively, the nucleation rate = l/Tpg' mcan ) for traffic breakdown at 
the bottleneck associated with (31) is 

g (B) _ ggU Fs(edAc) 1 / 2 / 8iV d A^ 



Note that <? determ) FS , N d , and £ d depend g on . Therefore, the mean time delay 
Tpg' mean ' ) (31) and the nucleation rate G^J ( 32) are functions of g sum and q on . 

If the flow rate q on decreases continuously up to a small enough value (how- 
ever, we assume that q on > 0, i.e., the deterministic cluster still exists at the 
bottleneck), then the values £ d , iV d , and q^ eim FS in (30), (32) and, therefore, 

( B) 

the nucleation rate G FS ( 32) do not decrease proportionally to this decrease in 
q OQ . In contrast, in this limit case £ d -> £ d , Um , N d -> iV di Um , and ^term, fs -» 
?determ, iim> where Hm, N d) Hm, and gJLm, lim are constants. Taking into ac- 
count that in this case in A c (30) the flow rate ? d et er m, fs ~ 9determ, Hm = const, 
we can see that at small enough values of q on the nucleation rate for traffic 
breakdown (32) depends on the total flow rate g sum only. In other words, in 
this limit case at a given q sum within the flow rate range (17) the nucleation 
rate for traffic breakdown at the bottleneck (32) tends to a finite constant 
value, which is greater than zero (see Sect. 4.1). 

When q on = 0, the road can be considered homogeneous one without bottle- 
necks. Then there is no deterministic perturbation (deterministic cluster) at 
the bottleneck and, therefore, the nucleation model and results of this article 
cannot be applied. In three-phase traffic theory, the breakdown phenomenon 
can also occur in this case. However, at the same conditions, in particlular, the 
same flow rates downstream of the bottleneck and on a homogeneous road, 
the nucleation rate for traffic breakdown on the homogeneous road is consid- 
erably smaller than at the bottleneck [33,3]. This is associated with empirical 
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results in which the breakdown phenomenon has also been observed away from 
bottlenecks, however, this traffic breakdown is much more rare than at an on- 
ramp bottleneck [3]. A consideration of a nucleation model of the breakdown 
phenomenon for a homogeneous road is beyond the scope of this article. 

As usual for each first-order phase transition observed in many other systems 
in natural science [38], the nucleation rate for traffic breakdown (32) is an 
exponentional function of A c (30). For traffic flow, in accordance with (32) 
and (30) the exponential growth of the nucleation rate with A c (30) is very 
sensible to the critical value for the deterministic breakdown phenomenon 
?determ fs- This emphasizes the important impact of the deterministic cluster, 
which occurs at the bottleneck at q on > 0, on the nucleation rate for traffic 
breakdown (32) at a given total flow rate g SU m- 



4 Discussion 



4-1 Numerical simulations of general results of nucleation model for traffic 
breakdown 



Let us compare general results of the nucleation model presented in Sect. 3 
with the diagram of congested patterns at an on-ramp bottleneck postu- 
lated in [30] and found in numerical simulations of microscopic traffic flow 
models [31,32], as well as with a microscopic theory of the breakdown phe- 
nomenon [33]. To reach this goal, we consider an example of the function W- 
(14) 



W -(N) = N 



i 

1 + (iV/iVo) 4 + . 



(33) 



where a, b, and N are functions of q OIl : a(q on ) = 1.32qo(q on )/N (q on ) 1/h, 
g (g on ) = 2700 + 370(1 + g on /300)- 1 vehicles/h, b(q on ) = SS + lofl + go^SO)- 1 
1/h, N (q on ) = 25 — 6.5(1 + g on /300) _1 vehicles; the unit of q on is vehicles/h. 

The analytical function (33) allows us to perform a numerical analysis of 
the mean time delay (20) and the associated nucleation rate (24) for the 
breakdown phenomenon (F— >S transition). For the analysis of (20) and (24), 
only branches N( dctcim ) and N c (Fig. 3) of the function (33) associated with 
the deterministic and critical clusters within which N < N s are relevant. This 
is because the maximum possible value of N = N 2 in the potential barrier A$ 
(21) that determines the nucleation rate (24) is equal to N s . 

However, for a qualitative illustration of a possible synchronized flow state re- 
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N [vehicles] N [vehicles] 



Fig. 5. N-shaped function <Zj q ° (N) (33) (a) for q on = 100 vehicles/h (curve 1) and 
q on = 600 vehicles/h (curve 2), and the associated potential (22) (b) as functions 
of the vehicle number N for q on = 100 vehicles/h and for three different total flow 
rates g sum : 2070 (curve 3), 2200 (curve 1), 2400 vehicles/h (curve 2). 

suiting from an F— >S transition, in (33) the branch for the synchronized flow 
state is added, in which the detachment rate W- (N) increases with N 3 . This 
branch corresponds to N > N s . Respectively, this branch of the detachment 
rate w-(N) has no influence on the analysis of (20) and (24). For this reason, a 
simple mathematical approximation (33) of the latter branch oiw-(N) is cho- 
sen, in which the detachment rate W-(N) exhibits formally unlimited growth 
with N. As mentioned, this does not impact on results discussed below. More- 
over, as follows from (A. 4) (see Appendix) probability of cluster emergence, 
which size N is large, is negligible. 

A numerical study shows that the potential $ exhibits qualitatively the same 
behavior at different total flow rates g sum (Fig. 5 (b)) as those in Fig. 4. The 
potential barrier A$ in (20), (24) (Fig. 6 (a)) and the associated critical ve- 
hicle number difference AN C (18) are decreasing functions of the total flow 
rate g sum ; at a given g sum they can also be decreasing functions of q on (Fig. 6 

3 For more detail explanation of the approximation (33), note that as in the general 
model (Fig. 3), the function (33) is a N-shape flow-vehicle- number characteristic 
(Fig. 5 (a)). This N-shape is chosen to satisfy those results of a microscopic three- 
phase traffic theory [3] in which a Z-shaped speed-flow characteristic for an F^S 
transition has been found (Fig. 1 (b, d)). The branch Vsyn, aver on this characteristic 
(Fig. 1 (d)) as well as the associated branch _/v( syn ) on the N-shape flow-vehicle- 
number characteristic (Fig. 4 (a)) are associated with a synchronized flow state, 
which results from the F— >S transition. The greater the density, i.e., the vehicle 
number N within the synchronized flow state, the greater the flow rate <Zdown 
(Fig. 1 (c, e)), i.e., the detachment rate w-(N) (14). In real traffic flow, the growth 
of W-(N) with N has obviously a limit. This limit is related to spontaneous moving 
jam emergence in synchronized flow of lower speed and greater density (i.e., greater 
N). In this case, an SP transforms into an GP, which consists of two traffic phases, 
synchronized flow and wide moving jams [3]. However, these effects are beyond the 
scope of this article. As mentioned above, the simple mathematical approximation 
(33) of the branch of i«_ (N) for synchronized flow states that are associated with 
N > N s can be used, because at N > N s the function w~{N) has no impact on 
results presented in the article. 
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Fig. 6. Potential barrier A$ (21) (a), critical vehicle number difference AiV c (18) 
(b), mean time delay for traffic breakdown Tpg' mean ) (20) (c, d), and nucleation rate 

for traffic breakdown G FS (24) (e, f ) as functions of the total flow rate q sum for three 
different flow rates q on : 100 (curves 1), 300 (curves 2), 800 (curves 3) vehicles/h. 



(b)). For these reasons, the total flow rate dependences of the mean time de- 
lay Tpg' mean ' ) (20) (Fig. 6 (c, d)) and of the associated nucleation rate for 
traffic breakdown at the bottleneck (Fig. 6 (e, f)) exhibit qualitative features 
observed in traffic flow at on-ramp bottlenecks [8,9] and found in a micro- 
scopic three-phase traffic theory [32,33]. This confirms that the breakdown 
phenomenon at the bottleneck is a first-order F— >-S transition [3] . In all these 
curves, the total flow rate g sum is smaller than the critical flow rate for the 
deterministic traffic breakdown ^eterm fs (4)- This means that A c (30) is not 
equal zero for all results in Fig. 6, i.e., traffic breakdown occurs due to a ran- 
dom density increase within an initial deterministic cluster at the bottleneck. 
The total flow rate dependences of the nucleation rate for traffic breakdown 
(24) calculated at different flow rates q on exhibit features of three-phase 
traffic theory in which the breakdown phenomenon can also occur at small 
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Fig. 7. Characteristics of the nucleation model: (a, b) - Boundaries of constant values 
of the nucleation rate G FS of traffic breakdown (curves 1-4), the critical boundary 



,-»(B) 

determ, S 



for deterministic traffic breakdown (curves F 



boundary 



fry\ 

curve F\ ) as functions of the flow rates q on and q m . (c, d) 



(B) 

determ, S^' 



and the threshold 



Dependencies of the flow rate (35) (curves 1-4), the critical flow rate detenu i s 

(4) for deterministic traffic breakdown (curves <^ eterm fs)> ana - * ne threshold flow 

rate g th (19) (curve g th ) as functions of q on . Curves 1-4 are related to different 

given values C for the nucleation rate of traffic breakdown G Fg in (34): 1/3.5 (curves 
1), 0.2 (curves 2), 0.1 (curves 3), 1/60 (curves 4) min . The nucleation model 
cannot be applied for q OQ = 0, therefore, the points in all figures in the vicinity of 
Qon = show only the tendency of the boundaries in (a, b) and the flow rates in (c, 
d) for the limiting case of small values q OQ in which, however, the on-ramp inflow 
rate q on > 0, specifically, it is assumed that the deterministic cluster still exists at 
the bottleneck. 



-(B) 



values q on . 

The critical boundary F$ £ (Fig. 7 (a, b)) in the diagram of congested patterns 
at the bottleneck (flow-flow plane with the coordinates (q on , q- m )) is associated 
with the cases in which the nucleation rate for traffic breakdown is a given 
value (. Therefore, the boundary F§ B ^ satisfies the condition 



G 



FS ((7sum> <7ony 



C, C = const, 



(34) 



i.e., at the boundary F$ £ the flow rate 



^^on) (35) 
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depends on q on . This boundary is qualitatively similar with the critical bound- 
ary Fg in the diagram at which the probability for traffic breakdown for a 
given time T Q b for observing traffic flow is 1 [3]. In the diagram, there is also 
the threshold boundary F th (curve F th in Fig. 7 (b)) at which the condition 
(19) is satisfied. The threshold boundary also exhibits the same qualitative 
features as those found in simulation of phase transitions and spatiotemporal 
congested patterns in a microscopic three-phase traffic theory [33,3]. In par- 
ticular, in the limiting case of small values q on (but q on > 0, i.e., it is assumed 
that the deterministic cluster still exists at the bottleneck) the flow rate 
reaches the maximum (limit) value <7q lim at a given nucleation rate for traffic 
breakdown ( (34). 

The greater the nucleation rate for traffic breakdown, the greater lim should 
be. However, the increase in nucleation rate for traffic breakdown has a limit 
associated with deterministic traffic breakdown occurrence: When ( in (34) in- 
creases, the boundary Fg B ^ for random traffic breakdown tends to the bound- 
ary F SLm, s (curves F d ( f t ) criIli s in Fig. 7 (a, b)) for deterministic traffic break- 
down in the diagram of congested patterns. At the boundary F^ erm s> the 
deterministic breakdown phenomenon occurs within the deterministic cluster 
even if no random vehicle number increase within the deterministic cluster ap- 
pears at the bottleneck. When ( in (34) decreases, the boundary Fg £ tends to 

the threshold boundary (curve in Fig. 7 (b)). In accordance with a 
microscopic theory [33], in the nucleation model the flow rate (35) (curves 
1-4 in Fig. 7 (c, d)), the critical flow rate ^ctcrm fs f° r deterministic traffic 
breakdown (curves q^J crm F s in Fig. 7 (c, d)), as well as the threshold flow 
rate q{^} (19) (curve q[^ in Fig. 7 (d)) can be the smaller, the greater q OQ is. 

4-2 Comparison with earlier nucleation models for traffic breakdown 

In [13,14,15] nucleation models for traffic breakdown for a homogeneous circu- 
lar road have been developed (see Fig. 11 in the review [40]). However, rather 
than traffic breakdown (F— >S transition), in [13,14,15] a nucleation theory 
for wide moving jam emergence in an initially homogeneous free flow (F^J 
transition) has been derived. Indeed, in final results of this probabilistic the- 
ory the vehicle speed within the vehicle cluster is chosen to be zero and a 
fundamental diagram for traffic flow with the vehicle cluster derived in the 
probabilistic theory [15] is qualitatively the same (see Fig. 48 in [40]) as those 
first found in [10] in a macroscopic theory of free flow metastability associated 
with F— >J transition. This fundamental diagram is confirmed by empirical 
results associated with wide moving jam propagation (see Fig. 17 in [40]). 
However, even on a homogeneous road, traffic breakdown is governed by an 
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F^S transition rather than by an F^J transition [3]. Thus, the nucleation 
theory of [13,14,15,40] does not describe traffic breakdown on a homogeneous 
road. 

The F— >S transition that can occur away from freeway bottlenecks is a very 
rare event [3]. This is because a freeway bottleneck introduces a spatial non- 
homogeneity in free flow at the road. The average speed within this non- 
homogeneity, which is permanently localized in a neighborhood of the bot- 
tleneck location, is lower and the vehicle density is greater than on the road 
away from the bottleneck [3]. This explains why in empirical observations 
traffic breakdowns are mostly observed at bottlenecks [7,8,9,3]. 

The first nucleation model for traffic breakdown at an on-ramp bottleneck has 
been suggested by Kiihne, Mahnke et al. [41,40]. In this model (see Chap. 
17 in [40]), a hypothesis of three-phase traffic theory about the sequence of 
the F^S^J transitions that governs phase transitions at the bottleneck [3] 
have been taken into account. In addition, in accordance with this theory [3] 
a random vehicle cluster, whose occurrence can lead to an F^S transition, is 
localized at the bottleneck [41,40]. 

However, in this nucleation model a random vehicle precluster that emerges 
from fluctuations is necessary. This precluster, which consists of one vehicle 
(n = 1), should occur in an initial hypothetical unperturbed free flow at 
bottleneck in which no vehicle cluster exists before. The precluster, which 
can be associated with a random decrease in speed of one of the vehicle in a 
neighborhood of the bottleneck, foregoes subsequent vehicle cluster evolution 
towards a critical cluster (critical nuclei) for traffic breakdown [41,40]. The 
attachment rate of precluster formation w+(0) is equal to the flow rate to the 
on-ramp q on [41,40]: 

w+(0) = q on . (36) 

The attachment rate (36) of vehicle precluster formation does not depend 
on the flow rate q in upstream of the bottleneck. However, in real traffic flow 
both flow rates q on and q in exhibit random fluctuations. The formula (36) 
should be associated with the basic model assumption that at q on = no 
vehicle cluster can randomly appear, specifically no traffic breakdown is pos- 
sible [41,40]. Apparently the assumption (36) leads to the nucleation rate for 
an F-^S transition at the bottleneck that is proportional to the flow rate to 
the on-ramp [41,40]. 

Whereas for a homogeneous road the model assumption for the necessity of 
precluster formation [13,14,15] is physically justified, this is not the case for the 
nucleation model at the bottleneck introduced in [41,40]. To explain this, note 
that in contrast with the model of a homogeneous road, at q on > and q in > 
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an initial state of free flow at the bottleneck is non-homogeneous regardless of 
fluctuations. This means that even no fluctuations would occur in free flow at 
the bottleneck, nevertheless free flow is non-homogeneous in a neighborhood of 
the bottleneck [3]. This is because two different flows permanent merge within 
the merging region of the on-ramp - the flow with the rate q on and the flow 
with the rate q m . Vehicles merging from the on-ramp onto the main road force 
vehicles on the main road to slow down. In turn, these slower moving vehicles 
on the main road force vehicles merging from the on-ramp onto the main road 
to decrease the speed too. Thus, the speed is lower and the density is greater 
at the bottleneck, i.e., a local cluster appears regardless of fluctuations. Thus, 
a permanent and motionless (deterministic) vehicle cluster in which speed is 
lower and the density is greater than away from bottleneck exists already on 
the road, even if there were no fluctuations in traffic flow. For this reason, the 
formula (36) [41,40] that assumes no vehicle cluster existence without random 
fluctuations in free flow at the bottleneck is in serious conflict with empirical 
results and results a microscopic three-phase traffic theory [3]. 

Moreover, the master equation of this model [41,40] searches the probability 
for a random vehicle cluster with n vehicles. In this master equation, a vehicle 
cluster exists (n > 0) only then, if the precluster has appeared. When there 
were no fluctuations at the bottleneck at all, then a vehicle cluster cannot ap- 
pear (n = 0) and no traffic breakdown is possible in the model. In contrast, in 
three-phase traffic theory deterministic traffic breakdown is possible regardless 
of fluctuations. This deterministic traffic breakdown occurs when the size of 
the deterministic cluster exceeds some critical value. This is the consequence 
of the non-homogeneity in free flow at the bottleneck mentioned above. 

Thus, in contrast with the above basic assumptions of Ref. [41,40], a nucleation 
model that can be adequate with empirical results should search probability 
p for random vehicle cluster evolution in which the cluster size N randomly 
changes due to fluctuations in a neighborhood of the deterministic cluster 
(Fig. 2). The size of this deterministic cluster ]\[( dctCTra ) does not depend on fluc- 
tuations in traffic flow. Random fluctuation either increases the speed within 
the cluster or decreases it. Consequently, the density and the cluster size N 
either decreases or increases. In the latter case, traffic breakdown occurs at a 
smaller flow rate g sum = q in + q on than the critical flow rate g sum = %et C rm fs 
associated with the deterministic traffic breakdown that occurs without any 
fluctuations at the bottleneck. 

These fundamental differences in the nucleation model of Ref. [41,40] and the 
model presented in this article can explain different results of these models. 
In [41,40], even if the on-ramp inflow rate q on is high, in an initial steady state 
of traffic flow there is no deterministic cluster at the bottleneck. As a result, 
there is no deterministic breakdown phenomenon of three-phase traffic theory 
in this model. Probably for this reason, in the model [41,40] the nucleation 
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rate for the breakdown phenomenon (generation rate of traffic breakdown 
critical nuclei) is proportional to the on-ramp inflow rate q on . As a result, if 
q on decreases below a small enough value (but q on > 0) and the total flow 
rate g sum increases (through an increase in q in ), a reasonable given nucleation 
rate for traffic breakdown at the bottleneck (the nucleation rate should be 
greater than m 1/20 min -1 , in accordance with empirical observations [8,9]) 
cannot be reached. This is true in the nucleation model of [41,40] even if the 
total flow rate is equal to a critical value associated with the critical nuclei for 
traffic breakdown consisting of one vehicle only (in [41,40] this critical value 
is denoted by q c2 ). 

In contrast, in our model the nucleation rate (generation rate of critical nuclei) 
for the breakdown phenomenon is not proportional to q on and, therefore, as 
mentioned in Sect. 3, for the limiting case of small values q on (however, we 
assume that q on > 0, specifically, the deterministic cluster still exists at the 
bottleneck) this generation rate of critical nuclei depends on the total flow rate 
<?sum only, i.e., this generation rate does not depend on q on . At a given q snm , an 
increase in q on can influence only on such characteristics of traffic breakdown 

fR^ (T$) 

as the critical flow rate <7 deter m fs anc ^ the threshold flow rate g th , as well as 
on congested traffic states at the bottleneck that result from the breakdown 
phenomenon. 

When g sum increases, the nucleation rate of traffic breakdown increases in 
the both models. However, in our model the nucleation rate cannot exceed 
the nucleation rate for traffic breakdown associated with the deterministic 
breakdown phenomenon. In contrast with assumptions of the nucleation model 
of Ref. [41,40], in our nucleation model the deterministic traffic breakdown 
occurs even without any random vehicle number increase within the initial 
steady state of free flow at the bottleneck. This is because if q on > 0, then in our 
model there is a deterministic vehicle cluster localized at the bottleneck, which 
exists permanent at the bottleneck due to the on-ramp inflow. In our model, 
random traffic breakdown nucleation can occurs through a random increase 
in vehicle number within this deterministic cluster. The mentioned qualitative 
differences in the nucleation model of Ref. [41,40] and our nucleation model 
are also responsible for different dependences of the generation rate of traffic 
breakdown on the total flow rate in these nucleation models. 



A Derivation of nucleation rate 



In order to derive formula (20), we use a general formula for the mean time 
delay T of escaping from the potential well for the master equation (6) [38] : 
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N 3 



n=Wi 



k=0 



(A.l) 



where p s (N) is a steady solution of (6), (7): 



N 

Ps (N)=p s (0)Yl—^- atiV>0. (A.2) 

n=l W -( n ) 

When Ni 1, i = 1,2 (more rigorous conditions are given below), the dis- 
tribution p s (N) has a sharp maximum at N — Ni, and the function p~ 1 (n) 
in (A.l) has a sharp maximum at n — N 2 . Then the formula (A.l) can be 
written as follows [38]: 



N 2 N 3 

T = K)" 1 J>(n) E Ps\n). (A.3) 

n=0 n=N 1 

Formula (A.2) can be written as 



p s (N) = p s (0) exp [-$(A)] at A > 0, (A.4) 

where the potential <3>(A) is given by (22). Substituting (A.4) into (A.3), we 
can find the exponentially large factor in (A.3) explicitly 



T = {w + Y l c lC2 exp [$(A 2 ) - $(A0], (A.5) 
where 

N 2 N 3 

Ci = E ex P ["A$ (1) W], c 2 = E ^p [A$( 2 )(n)], (A.6) 

n=0 11=^1 

A$ (i) (A) = $(A) - $(Ni), i = l,2. (A.7) 

The factors ci, c 2 can be estimated using the parabolic approximation of po- 
tential <£>(A) near the extremum points N — JVi, A 2 [39]. For instance, to 
find the factor Ci, we introduce a new variable y = N/Ni and approximate 
the sums in (22), (A.6) by integrals: 



A$« (A) « <jF> (y) = Ai / In ^ rfz, (A.8) 

N2/N1 

Ci«JVi J expf-^y)]*/. (A.9) 
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Using the series expansion 

4>W(y)=N 1 T h (y-l) 2 /2 + 0((y-l) 3 ) 



(A.10) 



near the point y — 1 , where rji = d\nw-/d\nN\ n=n 1 , we find c\ = ^j2irNi/r]i 



Similarly, c 2 = y ^N 2 /r]2, where r] 2 = —d\nw-/d\nN\ N=N2 . The substitution 
of Ci and c 2 into (A. 5) yields the formula (20). 

The parabolic approximation used for estimation of factor c\ holds only when 
we can neglect in integral (A. 9) third-order terms in the potential expan- 
sion (A. 10) [39]. The same is true for calculation of c 2 . The conditions of the 
parabolic approximation are 

Nirfi>il < = 1,2, (A.ll) 



where & = N 2 d 2 \nw_/dN 2 \ N=N ., i = 1,2. 

To derive the formula (28), we approximate the value A$ in (20) by integral 



A* = N 1 / ln W - {Niy) dy. (A.12) 

Under approximation of the function w_(N) in (A.12) near the maximum 
point N = N d by parabola, we get: 



W-(N) = ^term, FS 



1 , (g - iv d ) 2 

1 - U 



2Ni 



(A.13) 



where the formula (Zdeterm ps = w -(^d) is taken into account. The roots N = 
JVi and N = N 2 of equation g sum = ui-(N) given by formulae (14), (16) are 

N li2 = N d ^AN c /2, (A.14) 
where the critical value AiV c (18) is 

AN C = 2V2t;- 1/2 N d A 1 J 2 . (A. 15) 

Substituting (A.13)-(A.15) into (A.12), and calculating the derivatives w'_(Ni) 
in (23): w'^N,) » ±g£> FS £dAAy (27V d 2 ), i = 1,2, we get (28). 
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Under the condition that g sum is close to the critical point ^eiLm, fs> i- e -> when 
A c < 1, (A.16) 

from (28), we get the approximate formula (31). The latter is applicable under 
the condition (A.16) only if the conditions (A. 11) are still satisfied. Using the 
formula for derivatives w'_(Ni), % — 1,2 and that N i: % — 1,2 is close to N d , 
we can estimate the values & £d and ^ d AA r c /A r d , « = 1,2 in (A. 11). 
Thus, the condition (A. 11) reads ^AN^/N^ > 1. Taking into account formula 
(A. 15) for AN C , the condition (A. 11) can be written in terms of A c : 

N d A 3 J% 1/2 > 1. (A.17) 

This inequality together with (A.16) determine the range of A c in which the 
approximate formula (31) for the mean time delay Tpf' mcan - ) of an F— >S tran- 
sition at the bottleneck is valid. 
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